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ABSTRACT 

We investigate in details properties of stationary force-free magnetosphere of aligned 
rotator assuming the last closed field line lying in equatorial plane at large distances 
from pulsar. The pulsar equations is solved numerically using multigrid code with 
high numerical resolution, physical properties of the magnetosphere are obtained with 
high accuracy. We found a set of solutions with different sizes of the closed magnetic 
field line zone and verify the applicability of the force-free approximation. We discuss 
the role of electron-positron cascades in supporting of the force-free magnetosphere 
and argue that the closed field line zone should grow with time slower than the light 
cylinder. This yield the pulsar breaking index less than 3. It is shown, that models of 
aligned rotator magnetosphere with widely accepted configuration of magnetic field, 
like one considered in this paper, have serious difficulties. We discuss solutions of 
this problem and argue that in any case pulsar energy losses should evolve with time 
differently than predicted by the magnetodipolar formula. 

Key v^rords: starsimagnetic fields - pulsars:general - MHD 



1 INTRODUCTION 

Since the first works on pulsar magnetosphere, a stationary 
force-free magnetosphere of aligned rotator is considered as 
an underlying model for the real pulsar magnetosphere for 
more than 30 years. Despite its degenerated character (such 
"pulsar" even does not pulse), it is believed to reproduce 
qualitatively all main properties of the real pulsar magne- 
tosphere. For near-aligned pulsars it should give even an 
adequate detailed description. The structure of aligned rota- 
tor's force-free magnetosphere can be described by solution 
of a single scalar non-line a r PDE , t he so-called "pulsar equa- 
tion" , derived by l^^ch^l il973lJ l: [Scharlemann fc Wagoner! 
Jl973l) : IOkamotd JlQrST This is an equation for the flux of 
poloidal magnetic field. All other physical quantities describ- 
ing the magnetosphere are related to the fiux function <&, 
poloidal current J and angular velocity of magnetosphere's 
rotation Q, by algebraic relations. 

Analytical solution of this equation with non-zero 
poloidal current seems exist only f or the split-m onopole 
configuration of the mag netic field iMichellliggiTl and for 
a slightly perturbed split monop ole IBesHn et al.lll99d) . 
For dipole magnetic field analytical solutio n for the case 
of zero poloidal current has been found jMicheJ Il973al: 
iMestel fc Wana 19;^, but this solution is valid only inside 
the light cylinder (LC). There were several works dedicated 
to solution of linearised pulsar equation, where poloidal cur- 
rent and angular velocity were assumed to be proportional 



to magnetic flux function, what made the equation linear, 
but they did not lead to construction of a consistent model of 
aligned rotator magnetosphere (see e. g. Beskin ct al. IQSj 
lLvubarskii|[T990l : iBeskin fc Malvshkinl fl 998) . 

The flrst attempt to solve t h is equ ation numerically 
was made by IContopoulos et ahl 1^93), hereafter CKF. 
They have shown for the flrst time, that there exists a 
self-consistent solution with dipole magnetic fleld geome- 
try near the NS and magnetic fleld lines smoothly pass- 
ing through the light cylinder. In that work the position 
of the null point^ was fixed at the light cylinder and the 
question about applicability of the force-free approximation 
have been not investigated. Energy losses of the aligned ro- 
tator fo r CKF solution have been calculated by iGruzinovl 
J2005fl . iGoodwin et alJ J2004ll have studied this problem 
more deeply, namely they have searched for solution of the 
pulsar equation when the position of the null point is not 
fixed at the LC, but lies at different positions inside the light 
cylinder. For any position of the null point they obtained so- 
lutions, smoothly passing the light cylinder, but like CKF 
they have not studied physical properties of obtained solu- 
tions (e.g. energy losses, applicability of the force-free ap- 
proximations etc.). Their model, however, seems to be arti- 
ficial, because they assumed non-zero pressure in the closed 
field line zone, what implies continuous energy injection into 
the closed field line domain. Recently .Contopouloa L2005I) 
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^ the point where the last closed field line intersects the equatorial 
plane 
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addressed the case when the plasma rotation frequency in 
the open field line domain is different from the rotation fre- 
quency of the NS. It was shown, that there exist an unique 
solution of the pulsar equation for arbitrary plasma rotation 
frequency, although a rather simple case when the plasma 
rotation frequency is constant have been considered. Ap- 
plicability of the force-free approximation s in the magneto - 
sph ere of aligned rotator was considered in lTimokhinI (l2005l) 
and IContoDoulo j |2005). though in the latter work only for 
the null point located at the LC. 

Recently a different approach to the pulsar magneto- 
sphere modelling i s bein g developed b^Soitkovskv ( 20Q^, 
iKomissaro-JI ll2005l) . and iMcKinnevI 1I2OOQI . They perform 
time-depe ndent simulations of the pulsar magnetosphere. In 
Komissar ov (2005) aligned rota tor magnetosphere was mod- 
elled using full MHD code, in IMcKinnevI ll2006l) the same 
modelling was done with force-free code. The code of Ana- 
toly Spitkovsky allows performing of 3-D time-dependent 
simulation of the magnetosphere of inclined rotator by solv- 
ing equations of force- free MHD. In these simulations the 
existence of the stationary force-free magnetospheric config- 
uration was rigorously proved for the first time. Although 
this approach presents a big step towards the modelling of 
the real pulsar magnetosphere, in this paper we argue that 
it has serous limitation, namely the properties of cascades 
supplying particles in magnetosphere are not incorporated in 
these simulations. As it will be discussed later, cascades can 
set non-trivial boundary conditions on the current density 
in the magnetosphere. Its incorporation in time-dependent 
codes would require some efforts. 

In this work we investigate stationary problem solving 
the pulsar equation numerically with high numerical resolu- 
tion. We assume zero pressure in the closed field line region 
(cold plasma). As in all above mentioned works on numer- 
ical modelling of stationary aligned rotator magnetosphere 
we assume a topology with the current sheet in open field 
line domain flowing in the equatorial plane, i.e. configuration 
with Y null point - see Fig. Q This type of magnetosphere 
topology had became de facto the "standard model" , so we 
study it in details and analyse its properties regarding many 
aspects of electrodynamics. Smooth solutions are obtained 
for any position of the null point inside the light cylinder. 
High numerical resolution allows accurately incorporation 
of the return current flowing along the separatrix into nu- 
merical procedure. With the high resolution of used numer- 
ical method it was possible to calculate accurately physical 
properties of the solutions such as Goldreich- Julian charge 
density, magnetic fleld, energy losses and pointing flux dis- 
tribution etc., check applicability of the force- free approxi- 
mation and consider compatibility of the model with models 
of electron-positron cascades. 

Adjustment of the current density in the polar cap cas- 
cade zone of pulsar to the global magnetospheric structure 
was debated al ready in the first ten years after pulsar dis- 
covery (see e.g. lAronslll979|) . A concrete mechanism for the 
current density adjustment was proposed by Yu. Lyubarskij 
many years later, in 1992. At that time there was no self- 
consistent model of pulsar magnetosphere and detailed dis- 
cussion on this subject was difficult. Here we discuss the 
coupling between the polar cap cascade zone and the rest of 
the magnetosphere in the frame of the self-consistent model 
obtained in our simulations. We extend the picture proposed 
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Figure 1. Configuration of magnetic field in the magnetosptiere 
of aligned rotator witli Y null point — Y-configuration. After the 
null point xq the separatrix goes along the equatorial plane. The 
volume current / flowing in the open filed line zones [2] and [2'] 
closes somewhere beyond the light cylinder. There could be a 
volume return current along some open field lines, but the largest 
part of it fiows along the separatrix. 

by iLvubarskiil (Il992li addressing the evolution of the cur- 
rent adjustment mechanism with ageing of the pulsar. We 
also prove the necessity of such mechanism and discuss it 
in more details in the frame of the cascade model proposed 
bv lScharlemann et alJ lll97^ . We underline serious difficul- 
ties of model with Y null point regarding its compatibility 
with the Space Charge Limited Flow models of polar cap 
cascades and briefly discuss other possible magnetospheric 
configuration. 

The plan of the paper is the following. In Section |5| 
important properties of the pulsar equation are discussed. 
Model used in the current work and numerical method are 
described in Section |21 Results of numerical simulations are 
presented in Section 2] In section |^ we discuss the role of 
polar cap cascades for the global structure of the magneto- 
sphere, consider in details properties of the force-free mag- 
netosphere with Y null point, and highlight problems of the 
"standard" model of aligned rotator magnetosphere. A dif- 
ferent topology of the magnetosphere, with X null point, is 
briefiy discussed at the end of the section. We summarise 
the most important results in Section |S] 



2 THE PULSAR EQUATION 
2.1 General equation 

Here we adopt the wide used assumption, that the entire 
magnetosphere of the neutron star (NS) is filled with plasma. 
In some works starve d magnetosphere configuration hav e 
been debated (see e.g. lSmith et alJl200ll: IPetri et alJl2002li . 
where there are several separated clouds of charged parti- 
cles near the NS and no particle outflow, however there 
are indications that such configuration is unstable against 
diocotron instability dSpitkovskv fc Aronsl l2002: SDitkovs^ 
|2004|) . Plasma in the magnetosphere has to be non-neutral 



On the force-free magneto sphere of aligned rotator 3 



in order to screen the longitudinal (directed along magnetic 
field lines) component of the electric field, induced by ro- 
tation of the NS. In presence of the longitudinal electric 
field charged particles would be accelerated and their radia- 
tion will lead to copious electron-posit ron pair produc tion in 
super-strong magnetic field of pulsar 

(IStHrrQ^llil^), 

what 

finally results in screening of the accelerating field. 

Charge density necessary for cancelling of the longitudi- 
nal elect ric field, the so-called Gol dreich- Julian (GJ) charge 
density jGoldreich fc JulianllT969l) . near the neutron star is 
given by 



PG.T 



n B 
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(1) 



where 17 is angular velocity on neutron star rotation, B is 
magnetic field and c is the speed of light. Assuming that NS 
has dipolar magnetic field, the ratio of the particle kinetic 
energy density in the magnetosphere to the energy density 
of the magnetic field at the distance r can be estimated as 
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where e and rUe are electron charge and mass correspond- 
ingly, 7?NS ~ neutron star radius, 7 - Lorentz factor of accel- 
erated particles, Bo magnetic field strength in Gauss near 
magnetic poles of the star and P - period of pulsar rotation 
in seconds. All these quantities are normalised to their typi- 
cal values in pulsars. This ratio is small, less than 1 per cent, 
in the region with the size ~ 10'^ radii of the neutron 
star. It could remain small even further, but here magnetic 
filed deviates substantially from the NS's dipole field due to 
currents flowing in the magnetosphere, and Ekin/es can be 
estimated only after a self-consistent solution for the magne- 
tosphere structure is found. So, in a large domain surround- 
ing the neutron star, we can use force-free approximation, 
when particle inertia is neglected, and equation of motion 
takes the form 



pE + -^\3xB] 



0. 



(3) 



Hence, electric field E is perpendicular to the magnetic field 
B. Charge density p and current density j in eq. @ can be 
found from the Maxwell equations (we consider stationary 
problem) 

V ■ E = 47rp , (4) 

VxB = —3, (5) 

c 

With help of these equations equation can be written as 

(V-i;)B + [VxB]xB = 0. (6) 

In the force- free electrodynamics (FFE)^ the only possible 
motion of charged particles across magnetic field lines is the 
drift in crossed electrical and magnetic fields with the veloc- 
ity 

ExB 



(7) 



^ hereafter we use this shorter name for force-free degenerate 

electrodynamics (see e.g. iKomissarov 2002: Blandford. 2002,. and 

references there) 



Obviously |C/d| must be less than c, or equivalently E must 
be less than B. Generally speaking, eq. @ can have solu- 
tions where |[/d| > c. The surface, where \U-d\ reaches c, 
is commonly referred as the Ught surface. Beyond the light 
surface, where \Ud\ > c, the force-free approximation can 
not be applied. FFE is not self-consistent, because particle 
dynamics is ignored. Hence, each solution of eq. ® should 
be always checked for applicability of the force free approx- 
imation. 

In the axisymmetric stationary case considered here 
magnetic field in cylindrical coordinates {uj, (j), Z) can be 
written as 



V* X erf 



47r / 



(8) 



where is the unit azimuthal, toroidal vector. In compo- 
nents 

(B^,B^,Bz) = (--az*, — -, -a„*). (9) 

Scalar function '5 is related to the magnetic flux <I>niag trough 
the circle with the centre at the point (0, Z) and radius tjj 
by <I>niag = {zu , Z) . So, lines of constant 'i' coincides 
with magnetic field lines. It could be easily verified, that in 
force-free case the scalar function I{vj,Z) is constant along 
magnetic field lines, i.e 



/(*). 



(10) 



/ is related to the total current J outflowing trough the 
above mentioned circle by J = 2nI{zu,Z). 

In the quasi-s tationary cas e the time derivative of B 
takes the form fsee lMestellll973h 



dB 
'dt 



= Vx([nxr]xB) . 



Substituting this into the Faraday's law 
1 , 



VxE 



-dtB, 



we get for the electric field 

. . X __ ^ , 1^ __ _ __ ^ _ 

E = xB - W = V*- Vy , 

c c 



(11) 



(12) 



(13) 



where V is the non-corotational (see below) part of elec- 
tric potential. The first term in 1131 is poloidal and only 
the second term could make a contribution to the toroidal 
component. In axisymmetric case d^V = and, hence, E 
is poloidal. In the force- free case E B, from this follows 
that E ■ (V'l'xe^) — 0. Consequently, E oc V* and we can 
write 

£ = -^V*, (14) 
or in components 

{E^,E4.,Ez) = (-^9^*, 0, -^az*) . (15) 

Substituting this expression together with eq. JHJ into the 
formula for the drift velocity Q, we get 



471 mp „ ^ 
(7 D ~ iiFZue^ B = Of xr — kB . 



(16) 



C £2 

So, the particle motions is composed from rotation with the 
angular velocity Qf and gliding along magnetic field lines. 
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Hence, fip is the angular velocity of magnetic field lines ro- 
tation. By substitution of eg. 1141 into the stationary Fara- 
day's law one find Vf2pxV4' = 0. This implies that Qf is 
constant along magnetic field lines; 



(17) 



Equation 1171 is the well known Ferraro isorotation law. 

Finally, substituting E and B from eqs. (|HJ, I14II into 
equation lO we get 



1 - 



V"* 



+ 



2 dl_ 



(18) 



This is Grad-Shafranov equation for the poloida l mag- 
netic fi el d, the so-called pulsar equation, derived by l^chell 
l^l973b^ : IScharlemann fc Waeoneil l|l97.'t l: lOkamotol l)l974) . 

This scalar PDE is of elliptical type. It is the poloidal part 
of the vector equation @. The toroidal part of eq. Q is 
simply the relation HlOjl . Pulsar equation has two integrals 
of motion - I and fip. If we know them, we can solve this 
equation for function 'I' and determine the poloidal magnetic 
field. Electric and magnetic fields and all other parameters 
of the force-free magnetosphere can be found, because they 
are connected to ^I*, J and fip by algebraic relations. In the 
frame of FEE / and fip are free parameters. They could be 
determined self-consistently in the full MHD, if also electro- 
magnetic cascades , setting bou ndary conditions, are taken 
into account (see lBeskinll200i'Jl . Nevertheless one can get 
useful results in the force-free approximation. Equation l|18|l 
has one singular surface, the so-called light cylinder (LC), 
where zu = c/f2p (\['(ci7, Z)). As it will be shown in the next 
subsection the difference between Qf and is small and the 
singular surface has shape close to a cylinder with the radius 
of Rhc = c/il. 

We normalise variables and Z to i?LC and introduce 
new dimensionless coordinates x = lu/R-^c and z = Z/Ri^c- 
We will consider the case of dipolar magnetic field on the 
NS. So, near the star the magnetic field is given by 



(^2 + ^2)3/2 



*0 



(x2+22)3/2 ' 



(19) 



where n — BoR'^s/2 is the magnetic moment of the NS 
and ^0 = fJ,/RLC- We normalise 'if to 'I'o and introduce 
dimensionless function tp = 4'/'I'o. Instead of poloidal cur- 
rent function I we introduce dimensionless function 5 = 
(47r/c)(_RLc/^o)^- Angular velocity of magnetic field line 
rotation is normalised to the angular velocity of the NS by 
the relation Qf{x,z) = l3{x,z)Q,. For these dimensionless 
functions the pulsar equation I18II takes the form 

2 2 /5^2;2 + 1 

At the light cylinder the coefficient by second derivatives 
goes to zero and the pulsar equation has the form 



Let us now discuss properties of functions fip and S. 



(21) 



2.2 Q,F 

From relations 11711 it follows that V is constant along a 
magnetic field line. Hence, we could rewrite eq. 1131 in the 
following form 



dV\ 



Comparing this expression with eq. p4^ we get 



fip 



+ 



dV 



(22) 



(23) 



If there were no potential difference between different mag- 
netic field lines and between them and the surface of the 
pulsar, Q,F were equal to fi. But, independently of NS 
surface properties, a potential difference along open mag- 
netic fi eld lines will be always build in polar cap region of 
pulsa r ij^id^miM^^^u^erlandl Il975l : IScharlemann et alJ 
Il978l: iMuslimov fc Tsvganlll992ll . This lead to formation of 
a particle acceleration zone, where force-free approximation 
is not valid and charged particles are accelerated by the 
longitudinal electric field. Electron-positron pairs produced 
in the strong magnetic field of pulsar by photons, emitted 
by accelerated particles, screen the accelerating field, and 
as pair-production rate grows very rapidly with the dis- 
tance, acceleration zone terminates in a rather thin layer 
called pair- formation front (PFF). Above PFF accelerat- 
ing field is screened and FFE can be applied. The size of 
the acceleration zone is small compared to the overall size 
of the magnetosphere, its height varies from ~ 100 m for 
young pulsa rs in model with no particle e scape from the 
NS surface iRuderman fc Sutherland! Il975l) to 1-2 stellar 
radii in models, where particles freely escape the star sur - 
face iSc harlemann et al.l[l973 : iMuslimov fc Tsvganlll992f) . 
Geometrically this small region could be neglected in the 
modelling of the global magnetospheric structure. The po- 
tential difference between NS surface and magnetic field lines 
should be taken into account by boundary conditions on V , 
which can be reformulated as boundary conditions on Q,f- 
Potential difference along a magnetic field line in the accel- 
eration zone is determined by the position of PFF, which 
depends on local geometry of magnetic field, close to the 
NS surface, and kinetic processes in the electron-positron 
cascade. 

By the order of magnitude the relative difference of ro- 
tation velocities of plasma and NS can be estimated as 



Q.~Q,F 



^ P Vrot 

2tt rpc 
2.28x10" 



11 



Bo 



(24) 



P^AV, 



where AV is potential difference between NS surface and 
PFF (in esu units) , rpc — R'^gQ./c is the size of the polar 
cap; Vrot = cAV/{Borpc) is the linear velocity of plasma 
rotation relative to the NS surface in the acceleration zone 
- see eq. (31) in Rudcrman & Sutherland (1975). 

In the model with no particle escape from 
the NS surface the potenti al difference is given by 
iRuderman fc Sutherlandlll975L equation (23)) 



AV ~ 5.24xlO-'P 



-1/7 ^ 



Pc 



lO^cm 



4/7 



f- 



Bo 



-1/7 



(25) 
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where pc is the curvature radius of magnetic field lines. The 
potential difference is measured in esu units. Substituting 
these expression into eq. I|24|l we get 



x9xij) + zS^iJ) = 



106 cm 



J V1012G/ 



(26) 



We see, that for relatively young pulsars, with periods P < 
0.3 s, this ratio is very small, ~ 1 per cent. Even if the field 
line curvature radius is of the order of ~ 10* cm, typical for 
dipole magnetic field, for P < 0.1 s this ratio is ~ 2 per 
cents. 

For the model where particles freely escape the 
NS su rface we use estimations from iHibschman fc AronsI 
||200l|) . The potential diff erence in the acceleration zone 



(27) 
(28) 



Hibschman fc AronjIioOll eqs. (17) and (18)) 



AV 



9.87xWP~ 



l.llxlO^V 



'(i012g)^' 

3/2 / Bp \ 

UOi^Gy' 



,10i2G. 

Here h is the height of PFF above the NS surface in units 
of i?NS • The above estimations for accelerating potential are 
for the cases when when h > rpc and h < Vpc correspond- 
ingly. The potential differences are in esu units. The heights 
of PFF position due to photons emitted by non-resonant in- 
verse Compton scattering (NIC), curvature radiation (CR) 
and resonant inverse Compton scattering (RIC) of acceler- 
ated particles are given by 



kmc 



'NIC 



hcR 



hmc 



0.40 P 



0.12P 



VlOi^G^ 

1/4 ^ Bo 



0.68P 



1012 G 



-1/2 



1/2 



19/12 / Bp \ 

U012gJ 



-5/6 
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Vioi2Gy 



1012 G, 

-7/3 



pl/2 
P 



2/3 



fp: 



(29) 



(30) 



(31) 



(32) 



.1012G. 

see lHibschman fc Aron!j(l200j) . eqs. (34), (32), (42) and (37) 
correspondingly. Label "c" correspond to the model where 
the NS surface is colder than the polar cap of the pulsar, 
heated by the return current. Te is the temperature of the 
polar cap in units of 10® K. The radius of curvature of mag- 
netic field lines is factor fp times the radius of curvature of 



a dipole field, i.e. fp 



Pc 



dip 



~i/2 pj(9 .2xl0'^cm). 
i200ir ). in most pul- 



According to Hi bschman fc AronsI I 
sar the PFF height is set by non-resonant ICS photons. In 
high voltage pulsar, ones with the shortest periods - mil- 
lisecond and youngest pulsar with P < 0.3 s, the PFF is set 
by curvature photons. In both of these cases the resulting 
height of the PFF is larger, than the size of the polar cap, 
h > Tpc. Resonant ICS is important only for high field pul- 



sars, with B > 1.2x10" G, 
this into account we get 



in this case h <^ 



\it) 

V n 



o.09P/,r-(^)"^ 
0.027 p^/Vp^/^r-(^) 
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(33) 
(34) 



0.023 



0.034 pi/ Vp 77^''^ (y 



Bp 
1012 G 
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2xlO"G 
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(36) 




"NS g^x|> = 



l|) =l))(Xo) 



Figure 2. Calculation domain and imposed boundary conditions. 
See text for explanation. 



The temperature of the polar cap T due to the heating by 
return particles is of the order of 10® K. If the NS temper- 
ature is higher than this value, than formula 13311 should be 
applied, in the opposite case - formula 13411 . The tempera- 
ture of the NS surface depend of neutron star cooling model, 
and for rather young pulsar it should be higher than 10® K. 
So, formula 13311 is applicable for young, hot pulsar, where 
it gives for 5Q,/i^ ~ 0.01. Hence, in the model with free par- 
ticle escape, the ratio SQ/^l is of the order of few per cents 
for the majority of pulsars. 

We see, that 1 — /3 is of order of few per cents for most 
pulsars in the model with free particle escape and for young 
pulsars in the model with no particles escape. We restrict 
ourself considering only such pulsars, where 1 — /3 is small. 
Then the last term in pulsar equation 13811 is small in com- 
parison with other terms and could be neglected. In the rest 
of the paper we assume 

fip^n. (37) 

This assumption simplifies the pulsar equation 120L which 
now has the form 



(a:'-i)(a..V + 9..v) + — 



-d^tP ~ss' ^0, 



(38) 



where S' = dS/dil^. Nonlinearity in this equation is now 
present only in the term with the poloidal current function 
S. 



2.3 Poloidal current S 

In contrast to ilp, being set by kinetic processes in the po- 
lar cap, S depends on the global structure of the magneto- 
sphere. Both inside and outside the light cylinder the pulsar 
equation 1381 1 is a regular non-linear PDE of elliptic type. 
At the light cylinder this equation under assumption 1371 
has the form 



a, 7/) ^-ss' . 



(39) 



If function S is known, condition 1391 can be considered as 
a Neumann type boundary condition at the light cylinder. If 
boundary conditions are set both inside and outside the LC, 
the equation should have an unique solution in both regions. 
Generally speaking, for arbitrary function S solutions of the 
pulsar equation inside and outside the LC will not match. 



6 A. N. Timokhin 



lim V / lim • (40) 

So, a smooth^ solution is possible only for a specific func- 
tion S and the problem of finding a solution of the pulsar 
equations becomes an eigenvalue problem for the function 
S. 

The position of the light cylinder is not known a priory. 
For Qf different from SI it has a rather complicated form 
and, even if f2F(\[') as a function of is given by a model of 
the polar cap cascade, the position of the LC as a function of 
X and z has to be found self-consistently together with the 
solution of the pulsar equation. However, as it was stressed 
above, for most pulsar the deviation of the LC from a cylin- 
der with the radius c/f2 is of the order of few per cents or 
less. Hence, a solution of equation 13811 should give a very 
good approximation to the real magnetosphere of aligned 
rotator. 

The other open question regarding the poloidal current 
term in the pulsar equation is the topology of the magne- 
tospher e. In works of the Lebedev Physical Institute group 
(see e.g. lBeskin et al.ll993l:lBeskin fc Malvshkin"l998^ a ge- 
ometry with X null point have been assumed, hereafter X- 
configuration, see Fig. llSf a"). In that case the pulsar equa- 
tion should be solved in 3 different domains, separated by 
the current sheets. The positions of the point A(and A') is a 
free parameter of such model. Setting the positions of these 
points and the point xq one fixes the boundaries and gets 
a well posed, although complicated, problem. Such topol- 
ogy of the aligned rotator magnetosphere was criticised by 
iLvubarskiil il99Ci) . because the only source of the magnetic 
field in the magnetosphere is the pulsar itself, and in this 
case it is not clear what would be the source of the mag- 
netic field in the outer domain. The most frequently consid- 
ered topology of the aligned rotator magnetosphere implies 
an Y-like null point, hereafter Y-configuration, see Fig. 
In this case the only available free parameter in the model 
is the position of the null point xq- Fixing position of this 
point we fix the whole geometry of the magnetosphere. So, 
we have an elliptic equation with boundary conditions set at 
all boundaries of the closed domains with known positions of 
the boundaries. We wish to emphasise here, that the choice 
of the magnetosphere's topology is an additional assump- 
tion in the frame of stationary problem. In the following we 
investigate in details the force-free magnetosphere of aligned 
rotator assuming topology with an Y-like neutral point. 



3 NUMERICAL MODEL 

We solve equation I38II in a rectangular domain, see Fig. |5| 
The boundary conditions are the following. On the rotation 
axis (z-axis) 

1p{0, z) =0, ZnS < 2 < 2max . (41) 

At the equatorial plane, in the closed field line zone 

d^ip{x, 0) =0, a;Ns < a; < a;o , (42) 

following from the symmetry of the system. In the open field 
line domain 

^ if solution is continuous its smoothness follows from eq. 1391 . 
because 55" is the same at both sides of the LC. 



■ip{x,Q) ^ il){xo,Q), a;o < 2: < Smax , (43) 

i.e. the separatrix lies in the equatorial plane. Close to the 
NS the magnetic field is assumed to be dipolar, so for x = 

Xns, ^ 2 ^ 2ns and ^ a; ^ a;Ns, z = Zns 

^ix,z)^r'{x,z)^j^^^^^. (44) 

Magnetic surfaces should be c ome r adial at large distance 
from the NS, see 'ingraham' f 19 731). On the other hand, 
in calculations of Contopoulos et al.l 1^9^, where pul- 
sar equation was solved in the unbounded domain, with 
boundary conditions at infinity implying finiteness of the 
total magnetic flux, magnetic surfaces became nearly ra- 
dial already at several sizes of the light cylinder. Rather 
different outer boundary conditions, with finite magnetic 
fiu x inside the light cyl i nder at infinity, have been used 
bv ISulkanen fc Lovelace^ f 199(f). However, time- dependent 
simulations of JK omissarov (2005); Soitkovskv (2005) and 
iMcKimievI i2006i) provide strong evidence for correctness of 
outer boundary conditions when magnetic surfaces at large 
distances from the NS are radial. So, at the outer boundaries 
of the calculation domain for Q < x ^ a^max,^ = Zmax and 

X — a^max , ^ 2: ^ 2inax 

xd^il^ + z = 0. (45) 

At the light cylinder two conditions should be satisfied; 
(i) the solution should be continuous, 

tIj(x--^1~ ,z)^ip{x^l'^,z), (46) 

and (ii) the condition 13911 . These conditions together 
provide smooth tran sition through the LC. Following 
iGoodwin et alJ (|200^ we expand function ip at the LC in 
Taylor series over x implying continuity condition 14611 . By 
substituting the resulting expansion into the pulsar equa- 
tion 1381 1 and retaining the terms up to the second order we 
get the following approximation to the pulsar equation at 
the LC 

4 d^:,^{l, z)+2 a..^(l, z) = 9, [SS'{1, z)] . (47) 

This equation is nothing more than a reformulation of the 
smoothness conditions 14611 . 13911 valid for the first and sec- 
ond order terms in Taylor series expansion of ip. As the 
numerical scheme we have used is of the second order, this 
approximation, as well as its discretization, has the same 
accuracy as the discretized equation in the rest of the nu- 
merical domain. In course of relaxation procedure we are 
trying to satisfy the conditions 14611 . 1391 . i.e. we solve equa- 
tion 14711 at the LC instead of the original equation 13811 . 
which is singular there. Equation 13911 is used for determi- 
nation of the poloidal current term SS'{ip) along the open 
field lines. 

In the closed field lines zone, tp > ■f/'iast = tp{xo, 0), there 
is no poloidal current, so 5*5' = 0. The return current needed 
to keep the system charge neutral fiows along the separa- 
trix. In the open field line domain by setting the boundary 
condition 1431 the presence of an infinite thin current sheet 
is already incorporated into the solution procedure. How- 
ever, when the separatrix goes above the equatorial plane 
we have to model the current sheet. We assume that the 
return current is flowing along the field lines corresponding 
to the magnetic surfaces [ipia.st , V'last + dxp] . The total return 
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current flowing in this region is calculated by integration of 
the term SS': 




SS' dij) 



(48) 



We model the poloidal current density distribution over Tp 
in the current sheet '(/'last ^ ip ^ ipisst + dtp by an even order 
polynomial function going to zero at the boundaries of the 
current sheet 



S'W = A 



+ 



~2 



if 



(49) 



where constant A is determined from the requirement 
jj/'iaat+d'/' ^^^-j _ Q g^^j ^ integer constant. The 

pulsar equation is then solved in the whole domain includ- 
ing the current sheet. Although the current sheet cannot 
be considered as a force-free domain, but doing so we cal- 
culate correct the influence of the current sheet on to the 
global magnetospheric structure, though the obtained val- 
ues of the physical parameters inside the current sheet are 
fake. 

We developed a multigrid numerical scheme for solution 
of equations I38II and 14711 . These equations have been dis- 
cretized using the 5-point Gauss-Seidel rule. The coarsest 
numerical grid was constructed in the way, that the light 
cylinder is at cell boundaries. Each subgrid was obtained 
by halving of the previous grid. Cell sizes in the region 
X < 1, z < 1 are smaller in order to accurate calculate the 
current along t he separatrix. We used F AS scheme with V- 
type cycles fsee iTrottenberg et al.ll200lll . The Gauss-Seidel 
scheme was used as both smoother and solver at the coarsest 
level. At each iteration step both in the solver and smoother 
the new value of the poloidal current term SS'{l,z) was 
calculated from the relation 13911 at each point of the LC. 
Then a piece-polynomial interpolation of SS' in the interval 
(0,'!/)iast) was constructed and the return current distribu- 
tion was calculated according to the formulae 14811 and 14911 . 
Then for each point (a;, z) in the calculation domain the cur- 
rent term was calculated as SS'{x,z) = SS' {tp{x, z)), and 
the new iteration was started. So, we solved the pulsar equa- 
tion in the whole domain avoiding a very time consuming 
matching of the sol utions inside and outside the light cylin- 
der as it wa s done b y ContoDoulos ct al. ( 1999); ContoDoulos 
l|2005^ and lGruzinov (2005 ), though in ContQaoulos (.2005i) 
this matching procedure have been accelerated. As a start- 
ing conflguration a dipolar magnetic field everywhere was 
used. We did not encounter any problems with the conver- 
gence of the scheme for any value of xq , but for xq very close 
to 1 the convergence rate becomes essentially slower. Typ- 
ical number of points along each directions we used in the 
calculations was 3000 — 6000. 

We performed calculations for different values of nu- 
merical parameters in order to proof the independence of 
the results on the domain sizes (xmax, Zmax), the "NS size" 
(xtis, Zfis), the width of the current sheet dtp and the form 
of the current distribution (parameter k), as well as on the 
iteration procedure stopping criteria and number of points 
along both directions. Changes in convergence criteria and 
decreasing of the cell size from ones used in the most of our 
calculations did not produce relative changes in solutions 
greater that 10"*. In Table values of ipia.Bt and energy 



Table 1. Properties of obtained solution with xo=0.7 and xq 
approaching the LC for different values of numerical parameters 



numerical parameters: 



results: 



dtp 2k 








W 


xo=0.7 


0.03 2 


(8,7) 


(0.0667, 0.056) 


1.717 


1.864 


0.03 4 


(8,7) 


(0.0667, 0.056) 


1.712 


1.853 


0.015 2 


(8,7) 


(0.0667, 0.056) 


1.697 


1.821 


0.03 2 


(8,7) 


(0.0333, 0.028) 


1.720 


1.870 


0.03 2 


(16,14) 


(0.0667, 0.056) 


1.717 


1.864 


xo=0.99 


0.08 2 


(16,14) 


(0.06, 0.06) 


1.255 


0.977 


a;o=0.99231 


0.04 2 


(5,5) 


(0.0462, 0.0525) 


1.230 


0.939 



losses of aligned rotator W (see next section), obtained in 
computations with different values of listed numerical pa- 
rameters, are shown for xq = 0.7 and xq approaching the 
light cylinder. One can see, that with an accuracy of the 
order of few per cents obtained solutions are independent 
on particular values of the numerical parameters. Solutions 
with other a;o's have similar behaviour. 



4 RESULTS OF CALCULATIONS 

Our choice of the boundary conditions at the NS, equa- 
tion 1441 1. corresponds to the case when the dipole magnetic 
moment of the star /i is parallel to the angular velocity vec- 
tor n, /i||n. In this case the GJ charge density in the polar 
cap of pulsar is negative and there are electrons, which ffow 
away from the polar cap. The poloidal current S in the open 
field line zone is negative (see definition of the poloidal cur- 
rent eq. l|Hl)- In the case of anti-aligned rotator, i.e. /x is 
antiparallel to £7, all signs of the physical quantities related 
to the charge and current should be reversed. 

Calculations have been performed for the following val- 
ues of a:o:0.15; 0.2; 0.3; 0.4; 0.5; 0.6; 0.7; 0.8; 0.9; 0.95; 0.99; 
0.992. An unique solution has been found for each of the 
above a;o's. Let us consider in details physical properties of 
the obtained solutions. 



4.1 Poloidal current 

The poloidal current density S, calculated from the for- 
mula 1481 . does not deviate for more than ~ 20 p er cents 
from the values given by the Michel's solution jMichell 
Il973d) 



S = -tp[2 



tp 



(50) 



see Fig|l] The smaller xq, the smaller this deviation. The 
structure of the magnetosphere depends strongly on the 
poloidal current distribution. In solutions with xo ^ 0.6 
there is a domain in the open filed line zone, where vol- 
ume return current flows. However, only a small part of the 
return current flows there, the main part flows inside the 
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Figure 3. Global structure of the magnetosphere for xq = 0.992 - top figures, xq = 0.7 - middle figures, xq = 0.2 - bottom figures. 
Magnetic flux surfaces are shown by thin solid lines, the labelled vertical lines are contours of the drift velocity, grey area is the domain 
where the GJ charge density is positive. Dashed line separates regions with direct (above the line) and return (below the line) volume 
currents. The separatrix is shown by the thick solid line. On the left figures almost the whole calculation domain is shown, on the right 
figures - the central part of the calculation domain. Distances along x-axis (horizontal) and z-axis (vertical) are measured in LC radius 
RhC- 



current sheet. The size of this domain gets smaller with de- 
creasing of xo, and for xq ~ 0.6 the return current flows only 
along the separatrix, see Fig. |3 Qualitatively this property 
of solutions could be explained as the following. At the LC 
the condition 13911 is satisfied, so if dx^p changes the sigh the 
same occurs with the current term SS' , and the poloidal cur- 
rent density changes the sign. Magnetic field lines close to 
the null point are bend to the equatorial plane, but at large 
distance they become radial. So, for xq close to 1 dxtp < for 
some field lines, and volume return current must fiow along 
them. When xq decreases, more an more magnetic field lines 
at the LC will be bend away from the equatorial plane until 
there will be no lines bend to the equator. For field lines 



bend from the the equatorial plane dxi^ > and there is no 
volume return current along them. 

A convenient representation of the current density in 
the closed field line zone could be given by the current den- 
sity distribution in the polar cap jpc. In our notations the 
current density in the polar cap of pulsar normalised to the 
Goldreich- Julian current density Jg.t = PgjC is given by (see 
appendix 1X1 ea. llA6ll ) 

ipc = l,?G,T| (51) 

where 6/6pc is the colatitude normalised to the colatitude 
of the polar cap boundary 9pc, it is connected to the func- 
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Figure 4. Poloidal current distribution in the open field line 
zone normalised to the poloidal current from the corresponding 
Michel's solution. Inside the current sheet, for i/'iaat ^ 'A V'last + 
dij), (not shown here) the poloidal current decreases to 



tion ip through the relation 9/9-pc ~ y^i^/ipuat- In Fig. |H1 
jpc is shown for several solutions with different a;o's. The 
current density never exceeds the corresponding GJ current 
density and goes to zero at the polar cap boundary. The 
latter property is the consequence of the assumed magne- 
tosphere's topology. Indeed, from the condition at the LC, 
eg. 13911 ■ the current density along a given magnetic surface is 
proportional to the partial derivative dxip at the LC, but in 
configurations with Y null point dxil^ = for tp — ipiast- The 
deviation of the current density jpc from the GJ current den- 
sity increases close to the polar cap boundaries with increas- 
ing of xo- For solutions with xq ^ 0.6 the current density 
jpc changes the sign at some point near the boundary. On 
the other hand, jpc never exceeds the corresponding Michel 
current density and approaches jMichci when xo decreases. 



4.2 Drift velocity and force-free approximation 

The drift velocity in our notations is given by 
\Ud I S7ti7 -Bpoi X 



ItD = 



B 



(52) 



1 + 



Bpoi is the poloidal component of magnetic field. The light 
surface, i.e. the surface where the force-free approximation 
breaks down coincide with the surface, where ud = 1. We 
verified the applicability of the force-free approximations 
in each case. For most of the cases calculations have been 
performed in the domain with Xmax = 8, Zmax = 7, but 
for Xo — 0.2; 0.7; 0.992 we performed calculations also with 
Xmax = 16, Zmax = 14. lu all cases the light surface is located 
somewhere outside of these domains, see Fig.^ The drift ve- 
locity distribution for solutions with xo close to 1 even at 
large distances from the null point differs significantly from 
one in corresponding Michel's solution (the solution with 
the same ip\ast), where the drift velocity is the function of 
only i-coordinate. On the other hand, when xo decreases, 
md approaches the values from the corresponding Michel's 
solution. 



4.3 Charge distribution in the magnetosphere 

Goldreich- Julian charge density in the magnetosphere in our 
notations is given by 



SS' 



■ dxtp 



Pgj = po 



Po 



4nRtc 



(53) 



Close to the rotation axis the GJ charge density is negative 
and with increasing of the colatitude it becomes positive. 
While for solutions with xo > 0.6 the domain of positively 
charged plasma extends to infinity, for solutions with smaller 
xo's it becomes finite (cf. plots for xo = 0.2 with other plots 
in Fig. EJ. The reason for this is the following. At large 
distance from the light cylinder magnetic field lines becomes 
radial, so dxip is always greater than 0. Hence, there only the 
term SS' is responsible for changing of the charge density 
sign. However SS' for xq ^ 0.6 never changes the sign, see 
left plots in Fig. |21 For the same reason the volume return 
current always flows trough the positively charged domain. 
Close to the NS it passes trough the layer where charge 
density changes the sign, see right plots in Fig. |31 At this 
layer the so-called outer-gap cascade s hould develop (see e.g. 
ICheng et al.lll976l : iTakata et al1l2004l) . 

The force- free solution flxes not only the volume charge 
density, but also the charge density of the current sheet. 
As the electric field at opposite sites of the current sheet is 
different, the current sheet must have nonzero surface charge 
density. In Fig. |^ we plotted the linear charge density E 
of the current sheet as a function of distance I along the 
separatrix 



(54) 



where a is the charge density of the current sheet. E repre- 
sents the total charge of a volume co-moving with particles 
flowing along the separatrix with the constant speed, emit- 
ted at the same time (either at the NS or at "infinity"). 
E = const would imply a constant velocity flow of particles 
of one sign. However for each solution E is non-monotonic 
function with discontinuity in the null point. Such compli- 
cated dependence of E on / implies some non-trivial physics 
connected with particle creation in the current sheet, which 
is discussed in the next section. 

This complicated dependence of the current sheet 
charge density is easy to understand if one consider the 
so-called "matching conditio n" at the separatrix. As it was 
shown bv iLvubarskii il990l) '^. at the current sheet the fol- 
lowing condition for electric and magnetic field in closed (c) 
and open (o) field line domains should be satisfied 



J^c ~ J^c — t^o 



(55) 



This follows from integration of equation 1381 across the cur- 
rent sheet. In the closed field line zone there is no toroidal 
magnetic field. As it follows from eqs. 1141 and ||HJ the elec- 
tric field 

E = xBpoi . (56) 
Substituting this equation into equation 1551 we get 



B 



pol, c 



B 



pol, o 



(57) 



" see also lOkamotol 0^1), eq. (69) 
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Figure 5. Current density distribution in the polar cap of pul- 
function of the colatitude. jpc is normalised to the 
Goldreich- Julian current density |jgjI and the colatitude is mea- 
sured in units of the polar cap boundary colatitude 6pc. 



Figure 7. The total charge inside the sphere of the radius R 
centred at the NS: charge of the NS + charge in the magneto- 
sphere obtained by direct integration of pcj- Qtotal is normalised 
to 0.5 fi/ Ri^c- R is measured in units of -Rlc- 
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Figure 6. S - linear charge density of the current sheet (see 
text) as a function of the distance I along the current sheet. S 
is normalised to 0.5 fi/R^Q. I is measured in units of i?LC- The 
points marks the position of the corresponding null point. Note 
the jump in the charge density at these points. The dotted line 
corresponds to S = 0. 



From this and equation l|56^ follows that Ec > Eo and the 
charge density in the current sheet between closed and open 
field line domain, 



4-K 



Ec 



(58) 



is always negative. On the other hand, from the symmetry 
of the system - the electric field in regions 2 and 2' in Fig.d 
has different directions, - the charge density of the current 
sheet in the open field line zone 



4- 



(59) 



is always positive 

The total charge of the system, i.e. the charge of the NS, 
the magnetosphere and the current sheet together must be 



Figure 8. Energy losses of the aligned rotator as a function of 
xq in units of the corresponding magnetodipolar energy losses 

ll^mdl- 



zero. The boundary condition l|45|l implies that the total fiux 
of electric field through the sphere of a large radius is zero, 
hence the total charge of the system must be zero. In Fig. |7| 
the total charge inside the sphere centred at the coordinate 
origin is plotted as a function of its radius. The total charge 
of the system goes rapidly to zero at large distances from 
the NS. This plot could be also considered as an additional 
test of the numerical procedure, as the conservation of the 
total charge is not incorporated into the numerical scheme. 



4.4 Energy losses 

Energy losses of the aligned rotator in our notations are 
given by the formula 



W = |W„,di 



'Ala 



SdiP, 



(60) 
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Figure 9. Angular distribution of the energy flux dW/du) nor- 
malised to — 'i/'fj^at l^mdl/(47r), See eqs. I64i and 1651 . The distri- 
butions shown here are taken at ij = 4i?Lc ^i^d correspond to 
their asymptotic forms, see text. 



where |Wmd| is the absolute value of magnetodipolar energy 
losses, here defined as 



IT^mdl = 



4c3 



(61) 



see appendix rSl eqs. HBTII . IBSL In the obtained set of solu- 
tions W is function of xq. With decreasing of xo the amount 
of open magnetic field lines increases and, as the poloidal 
current dependence on ^ does not changes substantially, the 
energy losses of aligned rotator increases with decreasing of 
xq, see Fig. |S] Obtained dependence of energy losses W on 
the position of the null point xo could be surprisingly well 
fitted by a single power law 



W{xo) 



-0.94xo'-°'' \W^i\ 



(62) 



This formula is similar to the one obtained from analyti- 
cal estimations using Michel current distribution (see ap- 
pendix [HI eg.dEgl'l 
2 



W{xo) 



(63) 



The angular distribution of the energy flux (see ap- 
pendix [HI eq. (THlt 'l 



dW 
duj 



S- 



{z dxi> - X dzi)) 



(64) 



4% X 

In Fig.|U]this distribution is shown for several solutions with 
different xo- The Poynting flux distribution quickly reaches 
its asymptotic form at distance from the null point of the 
order of 1 — 2 iiLC- For example, in the case of xo = 0.99 
distributions taken at _R = 4 and 7? = 14 differs by no 
more than ~ 3 per cents. For configurations with smaller 
a;o's this deviation is even less. The smaller xo the close the 
angular energy flux distribution to the angular distribution 
in Michel's solution 
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Figure 10. Total energy of electromagnetic field in two difli'erent 
volumes of fixed sizes as a function of xq. H is normalised to the 
corresponding value of S(a;o = 0.15). 



dW \Wrni\ ,2 .2 

-T- = -. V'last sm ( 

duj 47r 



(65) 



because for small xo the solution at large distances is very 
close to the Michel's solution. In spite of recent works on 
modelling of jet-torus structure seen in Crab and other ple- 
rions (see|K omissarov fc Lvu barskv "200lf; Bog ovalov et alJ 



120051) . we note that magnetosphere configurations with 
larger xo would stronger support development of instabil- 
ities due to more asymmetric energy deployment into the 
plerion, providing more pronounced disk structure. 



4.5 Total energy of the magnetosphere 

The total energy of electric and magnetic fields in the mag- 
netosphere H = J{B'^ + E'^)/{8TT)dV would give information 
which configuration the system tries to achieve, the config- 
uration with the minimal possible energy. Obviously for the 
obtained solutions we could calculate the energy only in a 
finite domain. Another problem is very rapidly increase of 
magnetic field in the central parts, as r"''. As the magnetic 
fleld close to the NS is dipolar for each configuration, we 
calculate the total energy in a domain excluding the cen- 
tral parts. In order to verify the independence of the result 
on a particular domain we calculate the total energy in the 
magnetosphere in two different domains for each solution. 
These domains are defined as 0.2 ^ x ^ 5, 0.2 ^ z ^ 5 and 
0.075 ^ a; ^ 2.5, 0.075 ^ z 2.5 The results are plotted as 
a function of xo in Fig. 1101 The total energy of the magneto- 
sphere increases with decreasing of xo, so the magnetosphere 
will try to achieve the configuration with the maximal pos- 
sible Xq. 



4.6 Solution with xo —> 1 

The special case of — > 1 has been considered by sev- 
eral authors, because it was believed to be the real con- 
figuration of a pulsar magnetosphere 



Contopoulos et a n Il999l : lUzdenskvl lioo j ~l Gruzmovl 



JI^ibMskii 



199C 



200E 



Komissarovl 120051) . This case is peculiar in the sense, that 



magnetic field in the closed filed line zone diverges in the Y 
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Figure 11. Poloidal magnetic field strength in the equatorial 
plane as the function of x. Bp^i is normalised to ^/R^,-,. By the 
dashed line the theoretical prediction for _Bpoi(a;, 0) is shown for 
solution with xq = 0.992, ea. l66i . 

null point. Indeed, from equation 1571 it follows that near 
the null point, when xq ^ 1 



B 



\S\ 



pol 



(66) 



While the presence of the singularit y was noted 1 
iLvubarskiil and lUzdenskvl (l200d) . iGruzinovl dm 



firstly realised that such singularity is admitted, as it does 
not lead to the infinite energy of magnetic field in the re- 
gion surrounding the null point. In Fig. 1111 the strength of 
the poloidal magnetic field along the a;-axis is plotted for 
different solutions. By the dashed line the relation 11661 is 
shown. We see that when xo approaches the LC the mag- 
netic field inside the closed zone begins to grow close to the 
null point. This increase is more pronounced when the thick- 
ness of the current sheet decreases. Agreement between the 
curve for xq = 0.992, dtp — 0.4 and the dashed line is quite 
goo d. 

IGruzinovl ^2005^ solved an equation for the separatrix 
in the vicinity of the null point xq — 1 and have found 
that the angle at which separatrix intersects the equato- 
rial plane should be 77.3°. In our calculations we found this 
angle to be « 78° for xo = 0.992, dtp = OM and « 70° for 
xo = 0.99, dip = 0.08. So, our numerical solution shows good 
agreement with the analytical one. Energy losses found by 
IGruzinovl ll200? fl are 1.0 ± 0.1, what quite good agrees with 
va lues for W from Tab led Value of V'last = 1-23 calculated 
bv IContopoulo j (jioO^ coincide with ones from T able U and 
is close to 'i/'iast = 1.27 obtained bv IGruzinovl i2005l'l . al- 
though both of these results have been obtained with codes 
having worse numerical resolution than the code used in this 
work. 



5 DISCUSSION 

It seems naturally to assume, that force-free configurations 
are energetically preferably in comparison with configura- 
tions where there are geometrically large volumes with par- 



allel electric field^. Accepting this, we conclude that mag- 
netosphere of a pulsar should evolve through a set of force- 
free configurations. It does not necessary mean that for a 
relatively short transition time the system could not be es- 
sentially non-force-free, but rather that the most time the 
magnetosphere of an active pulsar is force-free. 



5.1 Polar cap cascades and force- free 
magnetosphere 

In a force free configuration the current density distribu- 
tion is not a free parameter, it is set by the structure of the 
magnetosphere, for example, by the value of xq in the case 
of Y-configuration. However, the current in the magneto- 
sphere of pulsar is supported by electron-positron cascades 
in the polar cap, i.e. the most of current carriers are pro- 
duced in the magnetosphere and are not supplied from ex- 
ternal sources. Independently of neutron star crust proper- 
ties, i.e. whether or not charged particles could be extracted 
from the surface, in polar cap of young pulsars electron- 
positron cascades are devel oped filling the magnetosphere 
of the star wit h particles ifRudemiM^^Su^erlan dl Il975l : 
IScharlemann et al. 1978; Muslim ov fc Tsveanl Il99^ . Also 
these particles are necessary in order to support MHD like 
structure of the magnetosphere. The current in the magneto- 
sphere fiows trough this cascade region, hence, the cascade, 
which properties depend on local magnetic field structure, 
has to adjust to the global properties of the magnetosphere 
too, namely to the current density flowing through it. We 
focus here on the case of stationary cascades. The hypothe- 
sis about stationarity of the polar cap cascades, when tem- 
poral variations of the accelerating electric field over the 
whole polar cap is much less than the accelerating field 
itself is widely adopted (e.g. IPaughertv fc Hardind Il982l : 
iRuderman fc Sutherlandlll975^ . We briefly address also the 
case of essentially non-stationary cascade iLevinson et alJ 

liooi). 

As it was shown bv ILvubarskiil (^g^), for current ad- 
justment in the stationary cascades a particle inflow from 
the magnetosphere into the cascade region is required. The 
typical current density, self-consistently supported by sta- 
tionary polar cap cascades, is close to Jg.t . For current densi- 
ties, both larger or smaller, than the Goldreich-Julian one, a 
particle inflow is necessary. The source of inflowing particles 
needed for current adjustment could be outer gap cascades, 
operatin g at the surface wh ere GJ charge density changes 
the sign ICheng et alj|T97^ . On the other hand, inflowing 
particles could be provided by the pulsar wind, where some 
outflowing particles could be reversed back to the NS due to 
momentum redistribution or due to small residual electric 
field arisen as the magnetosphere tries to support a force- 
free configuration. However, the zone where particles could 
fiow toward the NS is limited by the light cylinder (see Ap- 
pendixjnj. So, the source of inflowing particles must be in- 
side the LC. 

For iRuderman fc Sutherlandl (Il975l) cascades, when 
particles can not be extracted from the NS surface and are 
produced in the discharge zone, the adjustment mechanism 
works as follows. Inflow of positrons increases the current 

5 however see e.g. ISmith et all i200ll) : [pitri et al.l i2002ft 
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density, inflow of electrons decreases it. In the first case the 
inflowing positrons decrease charge density in the Pair For- 
mation Front (PPF) and more electrons is necessary to ad- 
just the charge density to the GJ value. This additional elec- 
trons together with inflowing positrons increase the current 
density. When there is an inflow of electrons, less primary 
electrons are necessary in order to support the GJ charge 
density at the PFF. Inflowing electrons are turned back at 
the PFF, and compensate the inflowing electric current. The 
outflowing current is only due to the primary electrons from 
the discharge zone, so the current density is less than Jgj- 

If particles could almost freely escape from the NS crust, 
the pulsar operates in the so-called Space Charge Limited 
Flow (SCLF) regime and the current density can not be 
essentially less then jgj. Indeed, the charge density in the 
discharge region, below PFF, is close to pcj and accelerating 
electr ic field forces charges to outflow with relativistic veloc - 
ities (|Scharlemann et alJllOTl iMushmov fc Tsvganlll992l) . 
For cascades operating in SCLF regime the mechanism of 
current adjustment works similarly for inflowing positrons. 
The particle inflow could increase the current density, but 
not decrease it. Only when the accelerating electric field is 
almost completely screened, the current density could be sig- 
nificantly less than Jg.i- However, in order to screen this ac- 
celerating field, charged particles infiowing from the magne- 
tosphere must penetrate practically up to the NS surface, i.e. 
they must have Lorentz factors comparable to the Lorentz 
factors of particles accelerated in the polar gap. In other 
words, somewhere in the magnetosphere inside the LC there 
should be zone(s) where particles are accelerated as effective 
as they would be accelerated in in the polar cap. Either the 
accelerating field there should be comparable to the one in 
the polar cap or the size of this zone would be essentially 
larger than some NS's radii. Both seems to be inappropriate. 

In both of these cases in order to support volume re- 
turn current, flowing in the direction opposite to jgj, the 
accelerated field in the polar cap discharge zone must be 
completely screened and the particles filling the magneto- 
sphere along magnetic field lines with return volume current 
must be produced somewhere in the magnetosphere. The 
accelerating electric field in the polar cap zone, being pro- 
portional to the magnetic field strength, is much stronger 
than any possible accelerating electric field far from the 
NS. Hence, the presence of the return volume current in 
the force-free magnetosphere seems to be incompatible with 
the force-free configurations of the magnetosphere, because 
the acceleration of particles to the required Lorentz factors 
with much weaker electric field requires large non-force-free 
domain(s) in the magnetosphere. The situation with non- 
stationary cascades is poor investigated, currently there is 
only one work dedicated t o detailed studies of s ignificantly 
non-stationary cascades - iLevinson et all ll2005 tl. However, 
we see no way how it would be impossible to support both 
particle production in the polar cap cascade and an average 
current having opposite dire ction to the d irection of accel- 
erating electric field, see also lAron^ il97!j) . 

In our consideration we assumed that the GJ charge 
density in the polar cap does not deviate sub stantially from 
its canonical value iGoldreich fc Julianlll969l) 



This is the case when the boundary of the polar can be 
considered as equipotential, i.e. having very high conduc- 
tivity. However if its conductivity is very low and surface 
charge density distribution at separatrix in the polar cap is 
different from the one in force-free solution, the GJ charge 
density can substantially deviate from values given by for- 
mula 16711 . In this case the characteristic current density 
flowing trough the cascade region would be different from 
the canonical value of --{yiBo) / (2ti) and, in principal, it 
could approach the values required by the global magne- 
tospheric structure, i.e. the problem of current adjustment 
could be solved by modifying Pgj instead of adjusting the 
deviation of j from Jg-i- Let us analyse this possibility. The 
largest part or the whole return current flows along the sep- 
aratrix. It could be electrons returning from region behind 
the light surface" of ions outflowing from the NS surface 
(see e.g.^uitkovskv & Arons 2004). If there are electrons in 
the current sheet close to the NS, then substantial devia- 
tion of electric fleld from the force-free value will give rise 
to electron-positron cascades producing enough particles to 
make separatrix near equipotential. Only ions, which much 
hardly emit photons capable to produce electron-positron 
pairs could support essentially non-equipotential polar cap 
boundary. However, as it was mentioned before, each partic- 
ularly force-free conflguration fixes the surface charge den- 
sity distribution along the current sheet everywhere where 
it is applicable. Independently on detailed structure of the 
polar cap zone the surface charge density along the separa- 
trix between closed and open field lines is negative, i.e. it 
must be enough electrons there, or the magnetosphere will 
be not force- free, see section and Fig.|S| Although in the 
discharge zone above the polar cap force-free approximation 
is not valid, and arguments of section 1^31 cannot be directly 
applied to the current sheet at the polar cap boundaries, 
electrons must be there for the following reason. The current 
sheet is a region where force-free approximation is broken, 
at least in some places, for example in the null point, where 
the surface charge density is discontiimous. As the return 
current flows in the current sheet, the parallel electric fleld 
will be directed from the NS, accelerating electrons in the 
current sheet toward the NS surface. Hence, in the current 
sheet at the polar cap boundary there are electrons too and 
this boundary will be approximately equipotential. Conse- 
quently, the GJ charge density in the polar cap should be 
close to the canonical value 16711 and in order to support a 
force-free configuration of the magnetosphere a current ad- 
justment mechanism is necessary. 

For current adjustment high particle density in the mag- 
netosphere is required. Indeed, only a small fraction of all 
particles could be turned back to the NS. There must be 
enough inflowing particles for adjusting of the current den- 
sity in the polar cap, i.e. its number density should be of the 
order of pcj/e- Hence, the particle number density in the 
magnetosphere must be 3> PG,i/e- However almost all par- 
ticles in the magnetosphere are produced in the polar cap 
and outer gap cascades, and a rather complicated coupling 
between cascade regions and pulsar magnetosphere arises. 
The weaker the cascades, the less particles are produced 



(67) 



" current sheet is not a force-free domain and considerations from 
Appendix |C] are not valid here 
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there, the smaller deviation from the GJ current density 
could be supported. Hence, when pulsar becomes older, the 
number of particles created in polar cap and outer gap cas- 
cades is smaller and the maximal deviation of the current 
density from j'gj will be smaller. If the magnetosphere re- 
mains force-free, its configuration must be changed in order 
to adjust to the new allowed current density. However, this 
new configuration would result in different energy losses of 
the pulsar, i.e. the ratio of the real losses to the losses given 
by the magnetodipolar formula will be different from the 
same ratio in previous configuration. So, generally speak- 
ing, the evolution of pulsar angular velocity derivative will 
not follow the power law oc — fi^, as it is predicted by the 
magnetodipolar formula. 

In the case of non- stationary cascades there are evi- 
dence that no particle infiow into the cascade region may be 
necessary in order to su pport current densiti es both larger 
and smaller than j'gj fsee lLevinson et al.l200.'Th . However, for 
creation of "wave-like " pattern of accelerating electric field 
iLevinson et all2005l) . necessary for support of small current 
densities together with reasonable pair creation rate, high 
pair density is required. With ageing of the pulsar the maxi- 
mal achieved electric field and pair density will decrease and 
shorten the range of allowed current densities. This would 
lead to the evolution of the magnetosphere similar to the 
case with stationary cascades. 

Arguments presented here are based on qualitative anal- 
ysis of the polar cap cascade properties. In order to make 
quantitative predictions a more detailed investigation of po- 
lar cap cascades is necessary regarding stationarity, ranges 
of current densities supported without particle infiow from 
the magnetosphere, and stability of the cascades in presence 
of particle infiow from the magnetosphere. 

5.2 Configurations with Y null point 

Let us analyse the behaviour of the magnetosphere of aligned 
pulsar under assumption that the null point is always of 
Y type. Here again we mean this in a time average sense, 
i.e. we neglect p o ssible non-stationa ry processes (see e.g. 
lKomissarovll2005l: IContopoulo3 l2005h in the current sheet 
operating on small scales (<C -Rlc), like building of small 
plasmoids. If non-stationary variations of the current sheet 
remains small, the stationary solution should adequately de- 
scribe the properties of magnetosphere. The total energy 
of the magnetosphere decreases with increasing of 2:0, see 
Fig. IIUI Apparently the system will try to achieve the con- 
figuration with the minimum possible energy, when xq = 1. 
However, when the restrictions set by the polar cap cascades 
are taken into account the picture becomes more compli- 
cated. 

In solutions with Y null point the current density in 
the magnetosphere close to the polar cap boundaries is 
always less than the GJ current density, it does not ex- 
ceed the Michel current density, see section 14.11 Fig. |^ 
The current adjustment mechanism could adjust the cur- 
rent density to the values less than ji'gj for stationary cas- 
cade model with no parti c le esc ape from the NS surface, in 
iRuderman fc Sutherland! lll975l) model. If pulsar operates 
in SCLF regime the current density can not be essentially 
less then Jgj- Hence, force- free solutions with Y null point 
are possible if the stationary polar cap cascade operates 




Figure 12. Current density distribution in the polar cap of pul- 
sar jpc as a function of the colatitude. Normalisation of physical 
quantities is the same as in Fig. |5] The colatitude ranges where 
the current density deviation from jgj could be supported by 
particles produced in outer gap cascades are indicated by dashed 
lines. The current density at colatitudes where jpc is shown by 
solid line should be supported by particles reversed inside the 
light cylinder. 

in Ruderman- Sutherland regime, or if the cascade is sig- 
nificantly non-stationary, the latter case however demands 
more detailed investigations. On the other hand, for solu- 
tions with xo > 0.6, the current density jpc close to the 
polar cap boundary has different sign than j'gj and such 
force-free configurations are probably never realised. 

As it was mentioned in section ISTTl the infiowing parti- 
cles could be produced either in the pulsar wind or in the 
outer gap cascades. The outer gap cascade could operate 
at the surface whe re GJ charge density changes the sign 
JCheng et al J 1197611 . Only relatively small amount of open 
field lines cross this surface. Hence, particles produced in 
the outer gap cascades could not adjust the current density 
along all magnetic field lines. On Fig. 1121 we plot current 
density in the polar cap of pulsar for several xo and indicate 
by the dashed line the colatitudes where particle infiow from 
the outer gap cascade would be possible. The critical colat- 
itude, where particle infiow from the outer gap cascade is 
still possible, corresponds to the field line with the smallest 
tl^ passing the surface of pc.i = inside the LC. At other 
colatitudes reversed particles from the pulsar wind (from 
inside the LC!) are necessary in order to adjust the current 
density. 

In Fig. 13 one can see that the deviation of the current 
density from jgj although remaining large, becomes smaller 
with decreasing of xo. So, if the magnetosphere remain force- 
free, with ageing of the pulsar, the configuration should 
change to the one with smaller current density deviation 
from jgj. Hence, if the force-free magnetosphere preserve 
its topology, with slow-down of the neutron star the size of 
the closed field line zone becomes smaller. Immediately con- 
sequence of this is the increasing of electromagnetic energy 
losses respectively to the corresponding "magnetodipolar" 
energy losses according to equation H62^ , see Fig. |H1 If at 
some time we approximate the dependence of xq on the an- 
gular velocity of NS rotation by the power law 
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xo cx n« , (68) 

where ^ is in reality a (complicate) function of pulsar age. 
5 > because xo decreases when pulsar became older. Sub- 
stituting it into the formula for pulsar energy losses (1621 we 
get 

Vf/ocfi", a = 4- 2.065 (69) 
and for pulsar braking index 

n = 5^ = a- 1 = 3- 2.065C, (70) 

i.e. the breaking index is always less than 3! 

Let us speculate that configurations with Y null point 
are energetically preferable over all possible solutions (force- 
free and non-force- free ones) and the polar cap cascade oper- 
ates in Ruderman-Sutherland regime. Then as long as parti- 
cles produced in the cascade regions will be able to support 
necessary current densities the pulsar magnetosphere should 
evolve with time as described above, and at each moment of 
time the configuration should be stable. Indeed, due to re- 
connection of open field lines in the equatorial current sheet 
the magnetosphere tries to achieve the energetically most 
preferably configuration, with xq = 1, but weaker cascades 
could not inject enough particles into the magnetosphere 
and support larger deviation of the current density from jg-j- 
So, 2:0 at each moment of time correspond to the configura- 
tion with current distribution having the maximal possible 
deviation from jgj • The polar cap cascade zone is the part of 
the whole system which do not allow clos ed field zone to hav e 
the maximal possible size. Conclusions of ' Spitkovskvl i2005l) . 
[Komissarov (2005) and McKiimov (2006) about instability 
of all configuration with xo < 1 is the result of assumption 
about possibility of arbitrary current density distribution in 
the pulsar magnetosphere. 

We should note also another peculiarity in Y- 
configuration - the jump in the surface charge density along 
the separatrix in the null point, where the charge den- 
sity changes the sign, see section 14.31 and Fig. |S| The re- 
turn current flows along the separatrix: electrons to the NS, 
ions/positrons from the NS. The surface charge density in 
the current sheet has different sign before and after the null 
point. We note that there is a jump in the charge density, 
not just a continuous changing of the charge density like it 
takes place across the surface where poj = 0. What happens 
in the null point, how such jump in the charge density could 
be supported when there is a continuous particle flow car- 
rying the return current? Do an electron-positron cascade 
operates here? Both the magnetic field and soft X-ray ra- 
diation of the NS are too weak here and electron-positron 
pair creation is suppressed. Anyway the problem requires 
additional investigations. 

5.3 Alternatives to force- free Y configurations 

In the force-free magnetosphere with Y-like null point the 
deviation of the current density from jgj is always large < 
I Jgj — jMicheil, especially close to the polar cap boundaries. 
For older pulsar with weak cascades it would be problem- 
atic to adjust the current flowing through the polar cap to 
the required value. On the other hand, even if force-free 
Y-configurations are energetically more preferably over all 



possible solution, the stationary polar cap cascade operat- 
ing in SCLF regime does not allow current density necessary 
to support an Y-configuration even in young pulsars. In that 
cases the magnetosphere could become non-force- free. 

The possible alternative to a magnetosphere, becoming 
non-force-free already at the distance of the order of 7?lc 
from the NS, would be a force- free magnetosphere with X- 
like null point. The current density deviation from the GJ 
current density in the magnetosphere with X-like null point 
would be less that in the Y-configuration for the following 
reason. In X-configuration with xq <1 condition (I39I I must 
be satisfied in points at the light cylinder above A and be- 
low A' , see Fig. 1131 In that points > 0, consequently 
SS'{ii < -(/'last) > and j(V' ^ i>ia.s\) < everywhere, nei- 
ther changing sign nor approaching zero. In principle, for xq 
not to close to the light cylinder and points A, A' not too 
close to the equatorial plane the deviation of the current den- 
sity from jG.T could be made rather small allowing even week 
cascades to support the current density, because in this case 
only small correction to the current density would be neces- 
sary. It is not clear now do force-free X-configurations exists 
and how they would look like (this work is in progress), 
however there is no clear physical reason forbidding such 
possibility. 

For X-configuration the jump in surface charge density 
of the current sheet in the null point could be avoided or at 
least reduced in magnitude. Considerations from section 
can be applied to the current sheet separating regions 2(2') 
and 3 in Fig. 1131 If directions of the poloidal magnetic field 
in these regions coincides'^, the charge density in this current 
sheet could be negative too. Indeed, the electric field in re- 
gions 2 {E2) and 3 (iSs) close to the separatrix in that case 
has the same direction. If £3 > E2, the charge density of 
the current sheen is negative and for such configuration the 
charge density at separatrix does not change the sign, but 
even if E2 > E3 the positive charge density would be less 
then in the case of equatorial current sheet, when the elec- 
tric field has to change the sign. If the directions of poloidal 
magnetic field in these regions are different, the same prob- 
lem with the current sheet as in Y-configuration arises for 
any values of the electric field inside the domain 3^ 

X-conflguration was criticised (iLvubarskiil ligQCTl be- 
cause in the force-free case there is no source of external 
magnetic field, which could fill the region 3 in Fig. 1131 
However, ones formed, this region may be supported by 
the global current system in the magn etosphere. Recon- 
nection in the equ atorial current sheet (lKomissarovll2005l : 
IContopouloll2005^ could lead to instantaneous formation of 
magnetic loops, which would grow in size and form some 
kind of X-configuration. These loops would try to merge 
with the closed field line zone or fly away, but current dis- 
tribution in the force-free magnetosphere will not support 
such configurations, so the resulting X-configuration could 
be stable. How would force-free solutions with X null point 
look like is not clear now, may be there could be solutions 
with many X null points, i.e. when there are several islands 



The current sheet in such configuration is necessary, because in 
the force-free magnetosphere charged particles can not flow across 
magnetic field lines, so after the null point they should flow along 
a very thin layer too. 
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Figure 13. Configurations of magnetic field in the magnetosphere of afigned rotator with X-nuII point. (a) - after the null point xq the 
separatrix goes away from the equatorial plane and never intersects it. There are 3 regions with open magnetic field lines: [2], [2'] and 
[3].(b) - after the null point xq the separatrix goes away from the equatorial plane but than intersects it (several) times in points xi, . . .. 
There are several regions with closed magnetic field lines: [1],[3] and [3']. 



of closed field lines along equatorial plane (see Fig. I13f b'l'l. 
but such complicated system may be unstable. On the other 
hand, the configuration shown in Fig. Il^f a)^ can not be 
entirely force-free. Indeed, there is no poloidal current in 
zone [3], as it has to have the same direction in both hemi- 
spheres, hence, the magnetic field there is purely poloidal. 
This implies that the plasma here does not rotate and there 
is no currents in the force-free case which could generate the 
magnetic field. However, we could speculate that in zone [3] 
at some (large) distance from the null point, where force- 
free approximation is broken, a system of currents is build 
up, which generate the magnetic field also in the force-free 
domain of zone [3], close to the null point. In this case the 
magnetosphere may be force- free at several sizes of the LC. 

Time-dependent simulations of aligned rotator magne- 
tosphere probably could cl ear what k i nd of c onfiguration 
is rea lised. Simulation s of Soitkovskv (2005); iKomissarovl 
^2005^ and lMcKinnevI ^2006.) do not incorporate restriction 
on the electric current in the magnetosphere due to polar cap 
cascade. In time-dependent codes the restriction set by the 
polar cap cascades should be formulated in form of bound- 
ary conditions on the poloidal current density. This could 
be done, for example, by introduction of an artificial Ohm's 
law along open field lines in the "polar cap" of pulsar, i.e. 
in regions close to the NS surface, like 

j = max((7||£, j^^J , (71) 

where j^^^ corresponds to the minimal possible current den- 
sity along the particular field line allowed by cascades. o-| | is 
conductivity along magnetic field lines, specific to each par- 
ticular code. In order to set these restriction the knowledge 
of polar cap cascade properties is necessary, what is another 
very complicated problem. 

The possibility of non-stationary magnetosphere also 



^ such configuration has been considered in several works (e.g. 
iBeskin et al.lil993t ISeskin fc Malvshkinlll998h 



can not be excluded at the current state of research. In or- 
der to decide between possible configurations (e.g. force-free 
with X or Y null points, non-force-free stationary magneto- 
sphere or significantly non-stationary configurations) more 
detailed studies of polar cap cascades and stability of the 
current sheet are necessary. 



6 CONCLUSIONS 

We have studied in details stationary configurations of the 
force-free magnetosphere of aligned rotator with Y null 
point. Assumption about Y-configuration of the magneto- 
sphere is very popular and this case had demanded care- 
ful investigations. We find a set of force-free solutions 
parametrised by the position of the neutral point xq- Re- 
sults presented in this work for xq = 1 agree very well 
with ones obtaine d by other authors iCon topoulo J l2005l: 
lGruzinovl2 005: Ko missarovll2005l: ISpitkovsk v 2005). We cal- 
culated physical characteristic of obtained solutions, and 
analysed properties of force-free magnetosphere with Y-like 
null point. For solutions with xo close to 1, we found that 
despite similarly distributed magnetic surfaces at large dis- 
tances from the light cylinder, they differs sub stantially from 
the split monopole solution of lMichel il973br) regarding dis- 
tribution of physical quantities (drift velocity, energy flux 
distribution, etc.). When the null point lies well inside the 
light cylinder, solutions approach the Michel one, the agree- 
ment being better for smaller values of xq's. 

We analysed the role, which cascades in the polar cap 
play in formation of the overall structure of the magneto- 
sphere. Although its properties depends mostly on the lo- 
cal physics in the polar cap of pulsar, this cascade region 
sets serious limitation on the current density in the whole 
magnetosphere. In some sense the non-trivial physics of the 
cascade plays a role of complicated boundary conditions for 
MHD equations describing structure of the magnetosphere - 
arbitrary current density are not allowed. Changes in bound- 
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ary conditions influence the whole solution. We argue than 
not all possible Y-configurations can be realised. Moreover, 
the restrictions set by some cascade models questions the 
existence of stationary force-free Y-configuration. To our 
opinion, there are two problems with the force-free mag- 
netosphere in Y-configuration: (i) current density strongly 
deviates from the GJ current density (ii) the charge density 
along the current sheet has discontinuity and changes the 
sign. These problems could be avoided in stationary force- 
free X-configurations. 

We argue that with ageing of pulsar and decreasing of 
the cascade power, the magnetosphere must involve with 
time, i.e. it should change to the conflgurations where the 
deviation of the current density from jc.i will be smaller. 
In the case of force-free Y-configuration the closed field line 
zone grows slower that the light cylinder during pulsar slow- 
down. This leads to decreasing of pulsar breaking index be- 
low the value 3, predicted by the magnetodipolar formula. 
This effect is present in aligned rotator because of the cur- 
rent adjustment in the polar cap of pulsar. Similar behaviour 
should be present in conflgurations with X-likc null point 
too. In analytical model of Bcskin & Malvshkin ( 1998) it 
was shown that the minimum energy of the magnetosphere 
is achieved when xo approaches the light cylinder. Although 
their model is an oversimplification of the real problem, this 
results should be qualitatively true. X-conflguration with 
Xo = 1 requires strong deviations of j from jc-i , so with age- 
ing of pulsar Xo should decrease, increasing number of the 
open fleld lines, leading to increasing of energy losses and de- 
creasing of pulsar breaking index. This should be true also 
for inclined rot ator, at least for not very large inclination 
angle. Recentlv lContopoulos fc Spitkovskv ( 200^ proposed 
another explanation for breaking index of pulsar being less 
than 3. In their model it is caused also by shrinking of the 
closed field line zone, but they assumed this is due to differ- 
ences in characteristic time, at which magnetosphere reaches 
the new configuration due to reconnection of new magnetic 
field lines, and the time of increasing of the light cylinder 
radius. However the reason of such slow reconnection in the 
current sheet is not clear. 

Obtained solutions could be used for comparison with 
observations. For example, magnetic fleld lines are differ- 
ently twisted in different solutions, this could be compared 
with magnetic field geometry inferred from pulsar polari- 
sation measurements fsee iDvk s fc Harding 2004, and refer- 
ences there). However, the real pulsar magnetosphere could 
be non-stationary and/or non- force- free, and this issue could 
be verified only within more detailed time-dependent ap- 
proach. On the other hand, in the latter case observational 
manifestations of pulsar will be signiflcantly different from 
models with force-free magnetosphere, where main emission 
comes from polar cap and from the outer gap zones. 
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APPENDIX A: CURRENT DENSITY IN THE 
POLAR CAP 

The poloidal current density is given by (see|Beski3|2003) 



J pol 



AI X 



dl_ 



'pol : 



(Al) 



the latter expression was obtained by taking into account 
relation dlj. Substituting for i?poi the expression for the 
dipole magnetic field at the NS surface Bpoi = -Bo(er cos6' + 
(l/2)eesin^) and expressing / and through normalised 
quantities we get for the poloidal current density in the polar 
cap of pulsar 



(A2) 



where | jc-i | is the absolute value of the Goldreich- Julian cur- 
rent density in the polar cap 



(A3) 



, RlC . 2 



(A4) 



For dipole magnetic field in the polar cap (see eq. 14411 ) 
Rns 

The colatitude of the polar cap boundary is 9pc = 1.45 x 
10~'^ P~^^^y/tpis,st, P is period of pulsar in seconds. So, 
9^ < <C 1 and the term with sin^ d in ea. llA2ll can be 
neglected. From equation IIA4II we have relation between the 
colatitude 9 in the polar cap and the corresponding magnetic 
flux function tp 

Ip = ^last I 77— I • (A5) 



Substituting this relation into equation 1A2II we get 

jpc = liG,T| i5'(('^] Vla=t). (A6) 



APPENDIX B: ENERGY LOSSES 

The poloidal component of Poynting flux in the magneto- 
sphere of aligned rotator is 

Ppoi = -^[i;xB]poi = -— /Bpoi, (Bl) 
47r c 



the latter expression was obtained with help of eqs. lIHt. 11411 . 
Expressing iJpoi through 4', eq. @, we have for the radial 
component of the Poynting fiux 

mp 1 



(B2) 



where r = ^/m^ + Z^ . Energy losses trough a soling angle 
duj are dW = —r^Pr dcj. The angular distribution of energy 
losses are given by 

dW _ Kh r 

duj 



C VJ 

Using normalised quantities we rewrite this equation as 



(B3) 



dW 

dui 



s 



(B4) 



47r X 

where |Wnid| is absolute value of magnetodipolar energy 
losses, here defined as 



4c3 



(B5) 



Energy losses of aligned rotator can be obtained by in- 
tegration of equation IB3II : 



W - 



dW 
duj 



duo - 



(B6) 



factor 2 appears because energy is carried away by the 
Poynting fiux from both hemispheres. Using normalised 
quantities introduced at the end of the section r2. II this for- 
mula can be rewritten as 



W^^c Sdi^=\W^i\ Sdij. (B7) 

^LC Jo Jo 

Analytical formula for estimation of aligned rotator en- 
ergy losses could be obtained on the following way. Poloidal 
current S for each of obtained solutions does not deviates 
much from the Michel's current function, ea. l|5()^ . Substitut- 
ing this function into equation lB7t we get 

W^~-|^Lt|M/md|. (B8) 

Dependence of i/iiast on xq we could estimate using mag- 
netic flux function of the dipolar fleld, ea. l|19|l . Substituting 
ip'^^'^{x(,) = into equation 1B8II we get 



(B9) 



APPENDIX C: THE SIZE OF THE REGION 
WHERE PARTICLE INFLOW IS POSSIBLE 

Charged particles in crossed electric and magnetic flelds drift 
with the velocity U-d given by eq. 11611 . which can be rewrit- 
ten as 

U-D — ^2 i^pol^'P ^ -B^Bpol) . (CI) 

In general, charged particle in the magnetosphere could have 
two velocity components: one perpendicular to the magnetic 
fleld line (the drift velocity) and an additional component 
along magnetic field line U\\, see Fig. ICll With increasing 
of the distance from the NS Ud increases, see Fig. H The 
velocity of the particle U can not exceed the speed of light. 
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hence, the paraUel component of particle velocity far from 
NS must be also smaller. Magnetic field lines far from NS 
are strongly twisted and at some distance the particle radial 
velocity component will be positive, i.e. directed from the 
NS, for any direction of the parallel component. At these 
distances particles can flow only from the NS. Let us show 
that the size of the domain, where particles could fiow to 
the NS, corresponds to the size of the light cylinder. 

Prom Fig. lCll it is evident that the radial component of 
the total particle velocity is 

Ur = Uj,r + U\\r- (C2) 

The maximal value of the velocity component parallel to the 
magnetic field is 

C/|f ~ Vc' - Ul . (C3) 

If the radial component of the magnetic field Br > Q then 
for the azimuthal component of magnetic field B^ — —\B^\. 
For radial components of particle velocities we have 



Unr = ^\B^\Br (C5) 

From this the maximal possible velocity component in radial 
direction is 




Particles could fiow to the NS only if (7™'"' < 0, and this is 
possible if 



B^^~c B~ ^^^> 



^ < -f , (C8) 

i.e. only inside the light cylinder. For Br < we get the same 
result. The same restriction, namely that particles could 



